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Section  I:  Project  Summary 


1.  Overview  of  Project 

This  project  is  performed  under  the  Office  of  Naval  Research  program  on  Basic  and  Applied  Research  in  Sea- 
Based  Aviation  (ONR  BAA12-SN-0028).  This  project  addresses  the  Sea  Based  Aviation  (SBA)  virtual  dynamic 
interface  (VDI)  research  topic  area  ‘Tast,  high-fidelity  physics-based  simulation  of  coupled  aerodynamics  of 
moving  ship  and  maneuvering  rotorcraft”.  The  work  is  a  collaborative  effort  between  Penn  State,  NAVAIR, 
and  Combustion  Research  and  Flow  Technology  (CRAFT  Tech).  This  document  presents  progress  at  Penn  State 
University. 

All  software  supporting  piloted  simulations  must  run  at  real  time  speeds  or  faster.  This  requirement  drives  the 
number  of  equations  that  can  be  solved  and  in  turn  the  fidelity  of  supporting  physics  based  models.  For  real-time 
aircraft  simulations,  all  aerodynamic  related  information  for  both  the  aircraft  and  the  environment  are 
incorporated  into  the  simulation  by  way  of  lookup  tables.  This  approach  decouples  the  aerodynamics  of  the 
aircraft  from  the  rest  of  its  external  environment.  For  example,  ship  airwake  are  calculated  using  CFD  solutions 
without  the  presence  of  the  helicopter  main  rotor.  The  gusts  from  the  turbulent  ship  airwake  are  then  re -played 
into  the  aircraft  aerodynamic  model  via  look-up  tables.  For  up  and  away  simulations,  this  approach  works  well. 
However,  when  an  aircraft  is  flying  very  close  to  another  body  (i.e.  a  ship  superstructure),  aerodynamic  coupling 
can  exist.  The  main  rotor  of  the  helicopter  distorts  the  flow  around  the  ship  possibly  resulting  significant 
differences  in  the  disturbance  on  the  helicopter.  In  such  cases  it  is  necessary  to  perform  simultaneous 
calculations  of  both  the  Navier-Stokes  equations  and  the  aircraft  equations  of  motion  in  order  to  achieve  a  high 
level  of  fidelity.  This  project  will  explore  novel  numerical  modeling  and  computer  hardware  approaches  with 
the  goal  of  real  time,  fully  coupled  CFD  for  virtual  dynamic  interface  modeling  &  simulation. 

Penn  State  is  supporting  the  project  through  integration  of  their  GENHEL-PSU  simulation  model  of  a  utility 
helicopter  with  CRAET  Tech’s  flow  solvers.  Penn  State  will  provide  their  piloted  simulation  facility  (the 
VLRCOE  rotorcraft  simulator)  for  preliminary  demonstrations  of  pilot-in-the-loop  simulations.  Einally,  Penn 
State  will  provide  support  for  a  final  demonstration  of  the  methods  on  the  NAVAIR  Manned  Plight  Simulator. 

2.  Activities  this  period 

During  the  period  of  this  report,  the  unsteady  flow  over  the  SPS2  ship  model  has  been  investigated  using  Craft 
Tech’s  flow  solver  CRUNCH  CPD.  An  unstructured  tetrahedral  grid,  generated  by  grid  generation  package 
Pointwise,  has  been  provided  by  CRAET  Tech  researchers.  The  unsteadiness  of  the  flow  was  demonstrated  by 
the  velocity  histories  with  respect  to  the  time.  Results  have  been  validated  with  experimental  data  presented  in 
[1]. 

CRUNCH-CPD  is  a  multi-physics  simulations  tool  for  analyzing  complex  flow  problems  developed  by  CRAET 
Tech.  It  has  three  different  modules  specialized  on  fluid  and  thermal  problems:  incompressible,  compressible 
and  thermal  [2].  In  the  scope  of  this  report,  CRUNCH-CPD  Incompressible  module  has  been  used  for  the 
unsteady  airwake  calculations  of  the  SPS2  Ship  model. 

A  2500x2500x1000  ft.  computational  domain,  shown  in  Pigure  1,  was  generated  by  Pointwise  mesh  generation 
software.  The  ship  body  and  sea  surface  were  designated  as  viscous  and  inviscid  wall  boundaries,  respectively. 
The  subsonic  inflow  boundary  applied  to  the  inflow  and  far-field  regions.  The  domain  was  chosen  as  large  as 
possible  to  avoid  the  potential  effects  of  grid  size  on  the  results.  The  full-scale  (455x45x55  ft.)  SPS2  (Simplified 
Prigate  Shape  2)  model  is  used.  Since  the  size  of  the  helicopter,  which  will  be  used  in  future  studies,  is  relatively 
much  smaller  than  the  ship  structure,  a  finer  mesh  topology  should  be  created  over  and  behind  the  ship  deck  to 


capture  the  detailed  flow  characteristics  of  ship  airwake  [1].  In  order  to  that,  the  unstructured  grid  was  clustered 
over  the  flight  deck  and  behind  the  superstructure,  which  can  be  seen  in  Figure  2. 


Figure  1  -  Computational  Domain 


Figure  2  -  Grid  structure  around  SFS2  Ship  model 

The  flow  characteristics  of  the  large  ships  are  highly  turbulent  and  unsteady.  The  superstructure  of  the  ship 
produce  large  separated  regions  and  vortex  flows  over  and  behind  the  ship.  Because  of  these  features  of  the  flow, 
the  airwake  of  the  ship  can  only  be  captured  with  unsteady  simulations.  However,  steady  state  simulations  can 
capture  the  main  tendency  of  the  flow  as  averaged  over  a  long  period  of  time.  And  these  results  can  be  used  as 
an  initial  condition  for  the  unsteady  CFD  simulations  to  speed  up  the  calculation  process. 

A  steady  state  simulation  with  a  At=0.001  sec  time  step  was  performed.  The  parallel  computing  was  carried  out 
using  COCOA4  computer  at  Penn  State  with  32  processors.  The  solution  convergence  was  determined  by 
monitoring  the  history  of  the  boundary  flow  residuals.  Roughly,  2000  iterations  were  needed  for  the  convergence. 

The  simulations  were  performed  as  unsteady  laminar  (MILES  LES)  which  has  been  shown  to  be  adequate  for 
airwake  simulations.  Laminar  flow  assumption  saves  from  time  and  computing  power  for  the  airwake 
simulations.  The  simulations  were  carried  out  with  a  free-stream  velocity  of  Voo=  12  m/s  (25  knots  or  42.2  ft/s) 
and  zero  side-slip  angle  (p  =  0  deg)  . 

Figure  3  and  4  shows  the  velocity  magnitude  distribution  over  the  domain  obtained  from  the  steady  state 


simulation.  Three  different  flow  separation  regions  can  be  observed  from  the  results.  The  results  show  that 
steady  state  simulations  can  capture  only  the  beginning  of  flow  separation  forming.  Also  it  can  be  observed  that 
the  flow  farther  from  the  ship  body  is  perfectly  steady  and  the  grid  domain  is  big  enough  for  the  airwake  region. 


Figure  3  -  Steady  State  CFD  Case  WODSpeed  =12.87  Figure  4  -  Steady  State  CFD  Case  WODSpeed  =  12.87 
m/s  (23  Knots)  WOD  Angle  =  0,  XY  plane  m/s  (23  Knots)  WOD  Angle  =  0,  X-Z  plane 

After  having  acceptable  results  with  steady  state  simulation,  the  simulation  has  been  restarted  using  the  steady 
state  results  as  an  initial  condition  for  the  unsteady  CFD  calculations.  1000  iterations  (outer  loop)  with  At=0.01 
time  step  have  been  carried  out  to  skip  transition  region  for  the  turbulent  flow  development.  After  that,  6000 
iteration  with  At  =  0.01  has  been  carried  out  to  capture  ship  airwake. 

Figure  5  shows  the  velocity  magnitude  distribution  over  and  behind  the  superstructure  within  a  time  period  of  10 
to  70  seconds  of  simulation.  Shedding  vortices  from  the  superstructure  can  be  observed  in  Figure  5.  Results 
show  that  the  periodicity  of  the  shedding  vortices  was  captured  successfully  by  the  flow  solver. 
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Figure  5  -  Velocity  magnitude  distribution  between  10^^  and  70**^  seconds  of  the  simulation 

The  numerical  results  have  been  compared  with  the  available  corresponding  time-averaged  measurement  data  on 
the  off-body  planes  in  the  flight  deck  and  the  superstructure  regions,  obtained  from  Ref.  [1].  The  experimental 
data  is  obtained  from  wind  tunnel  tests  using  a  1:8.5  SFS2  ship  model. 


Figure  6  shows  the  comparisons  of  the  time  averaged  streamwise  non-dimensional  velocity  distributions  over 
flight  deck  between  the  measured  experimental  data  and  the  computation.  On  Figure  6,  the  agreement  is  pretty 
good.  Figure  7-8  show  the  quantitative  comparison  of  streamwise  and  vertical  velocity  distributions  over  three 
different  planes  on  the  rear  ship-deck,  respectively.  The  plots  show  the  measured  data  at  17  ft  above  the  deck. 
The  agreement  is  pretty  good  for  the  streamwise  velocity  component,  but  there  is  a  small  difference  between 
calculated  and  measured  vertical  velocity  distributions.  This  error  might  be  a  result  of  laminar  flow  assumption 
on  numerical  calculations.  Also  the  experimental  data  belongs  to  a  wind  tunnel  tests  performed  with  a  1:8.5 
SFS2  model,  which  can  lead  to  scalability  problems.  Even  though,  the  sideslip  angle  of  the  flow  is  zero  and  the 
ship  geometry  is  symmetrical,  there  is  some  asymmetry  in  both  computed  and  measured  results,  which  is  also 
observed  by  Zhang  et  al  [1]  and  S.A.  Polsky  [3]. 


Figure  6  -  a)Time-averaged  streamwise  velocity  distributions  over  flight  deck,  CFD,  b)Time- 
averaged  stream  wise  velocity  distributions  over  flight  deck,  experiment,  [from  1] 


Time-averaged  streamwise  velocity  over  deck  @z=35.0ft 


Figure  7  -  Time  averaged  streamwise  velocities  over  deck  at  different  planes 


Figure  8  -  Time  averaged  vertical  velocities  over  deck  at  different  planes 


Figure  9  shows  the  comparison  of  time  averaged  streamwise  velocity  distributions  between  experimental  data 
and  calculated  data  over  the  super  structure  of  the  ship.  The  plots  shows  the  measured  data  at  7.5  ft  above  super 
structure  of  the  ship.  The  numerical  results  are  in  a  good  agreement  with  the  measured  data.  The  quantitative 
data  is  given  in  Figure  10  and  11.  There  are  some  offset  between  the  calculated  and  the  measured  data.  This 
offset  might  be  because  of  the  free-stream  velocity  difference  between  CFD  and  experiment.  However  the  trend 
is  correct  and  the  numerical  results  are  acceptable  for  the  current  phase  of  this  research. 


Figure  9  -  a)Time-averaged  streamwise  velocity  distributions  over  super  structure,  CFD,  b)Time-averaged  stream 
wise  velocity  distributions  over  super  structure,  experiment,  [from  1] 


Figure  10  -  a)  Time  averaged  streamwise  velocity,  longitudinal  plane  (x  positive  toward  ship  bow)  b)Time  averaged 

vertical  velocity,  longitudinal  plane 


Figure  11  -  a)  Time  averaged  vertical  velocity,  forward  plane  b)  Time  averaged  streamwise  velocity,  aft  and  forward 

plane 


3. 


Significance  of  Results 


The  results  obtained  using  CRUNCH  CFD  flow  solver  will  be  implemented  to  GENHEL-PSU  in  the 
next  stage  of  this  study. 

4.  Plans  and  upcoming  events  for  next  reporting  period 

•  A  subdomain  region  will  be  extracted  from  the  full  scale  model. 

•  These  flow  solutions  will  be  integrated  with  GENHEL-PSU  to  develop  a  set  of  baseline 
simulations  of  the  utility  helicopter  operating  in  a  ship  airwake  with  one-way  coupled  flow 
solutions.  This  will  provide  a  baseline  with  which  to  compare  the  fully  coupled  solutions. 

•  Begin  development  of  fully-coupled  simulations:  In  fully  coupled  solutions,  blade  position  and 
aero  loads  are  transmitted  to  the  CED  code,  the  CED  code  then  calculates  a  velocity  field 
(including  the  induced  velocities  from  the  aircraft  airloads)  and  sends  these  velocity  values  back 
to  the  helicopter  simulation  model.  The  subsequent  airloads  and  dynamics  of  the  helicopter  are 
then  affected  by  the  evolving  external  flow  field.  In  this  sense,  the  CED  solutions  serve  the 
purpose  of  not  only  the  ship  airwake  effects  but  of  the  induced  flow  field  generated  by  the 
helicopter  main  rotor  (and  possibly  other  components  of  the  aircraft).  Induced  flow  in  the  rotor 
is  usually  modeled  by  a  lower  order  model  in  flight  simulations  (e.g.  finite  state  inflow),  but 
these  modules  will  be  replaced  by  CED  in  the  coupled  solitions. 

•  Initial  coupled  solutions  will  not  involve  ship  flow  fields.  Coupled  simulations  will  be 
performed  with  the  helicopter  hovering  in  an  open  domain.  The  helicopter  will  be  trimmed  and 
perform  an  extended  hover  using  the  controller.  The  performance  and  trim  of  the  helicopter  will 
be  compared  to  those  predicted  by  the  simulation  model  without  coupled  CED.  We  expect  to 
begin  development  of  these  solutions  in  June,  with  results  expected  later  this  summer. 
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